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Positive systems analysis via integral linear constraints 

Sei Zhen Khong, Corentin Briat and Anders Rantzer 


Abstract —Closed-loop positivity of feedback interconnec¬ 
tions of positive monotone nonlinear systems is investigated. It 
is shown that an instantaneous gain condition on the open-loop 
systems which implies feedback well-posedness also guarantees 
feedback positivity. Furthermore, the notion of integral linear 
constraints (ILC) is utilised as a tool to characterise uucertaiuty 
in positive feedback systems. Robustness analysis of positive 
linear time-varying and nonlinear feedback systems is studied 
using ILC, paralleling the well-known results based on integral 
quadratic constraints. 

Index Terms —positive systems, closed-loop positivity, robust 
stability, integral linear constraints 

I. Introduction 

Monotone systems are a class of dynamical systems that 
preserve a closed order relation from the input space and 
initial state to their state space and output space [1]. They 
constitute one of the most important classes used in math¬ 
ematical biology and chemical modelling. Since biological 
and chemical models often treat variables such as popula¬ 
tion densities or concentrations of chemical mixtures that 
are intrinsically positive, they typically preserve positivity 
of solutions to differential equations and enjoy additional 
monotonicity or order-preserving properties. 

Positive monotone systems have gained an increasing 
attention over the last decade due to the fact that their 
order-preserving properties can often be exploited to simplify 
computations and control synthesis. For instance, stabilising 
output feedback controllers for positive systems are charac¬ 
terised using linear programming in [2] and extensions to 
input-output gain optimisation are considered in [3]. In [4], 
it is shown that the input-output gain of positive systems can 
be evaluated using a diagonal quadratic storage function and 
this is utilised for Hqo optimisation of decentralised con¬ 
trollers in terms of semidefinite programming. Rantzer [5] 
establishes that positive systems are amenable to distributed 
control whose complexity scales linearly with the number 
of interconnections. Robust stability of positive monotone 
systems subject to time-varying delays is studied in [6]. The 
class of differentially positive systems is introduced in [7] 
as an extension of the definition of monotone systems. 

Of the aforementioned works, [3] introduces the use of 
integral linear constraints (ILCs) to characterise uncertainty 
and derives linear-programming conditions to verify robust 
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stability of LTI systems based on a dissipativity approach. 
This paper further generalises this idea via an input-output 
approach to analyse robust stability of positive feedback 
systems, in a similar spirit to integral quadratic constraints 
(IQCs) based analysis [8]. IQCs are a well-known and useful 
tool for describing uncertain input-output dynamics. It unifies 
a range of methods in the literature, such as small-gain and 
passivity. This paper presents a first attempt to inherit these 
features via the more scalable form of ILCs for positive 
feedback systems. The main result generalises the static 
gain conditions on robust stability of positive feedback inter¬ 
connections of linear time-invariant (LTI) systems to linear 
time-varying (LTV) and nonlinear systems via integral linear 
constraints. Importantly, a sufficient condition guaranteeing 
positivity of the feedback interconnection of two positive 
nonlinear systems is also provided. 

The paper evolves according to the following stmcture. 
Notation and preliminary material are introduced in the next 
section. Section III establishes a sufficient condition under 
which positive feedback interconnections preserve positiv¬ 
ity of open-loop nonlinear systems. ILCs are reviewed in 
Section IV. Robust stability analysis via ILCs is performed 
in Section V. Section VI contains a couple of illustrative 
examples. Some conclusions are provided at the end. 

11. Notation and preliminaries 
A. Signals and systems 

Denote by L" the set of ]R”-valued Lebesgue integral 
functions: 

L” := : [0, oo) —> M" : ||n||i := J |v(f)| df < oo j-, 

where [ • | denotes the 1-norm. Let M” denote the nonnegative 
orthant of M". Define 

:= {w G Li : v{t) G M!; a.e.}, 

where “a.e.” is with respect to the Lebesgue measure on 
[0,oo). Define the truncation operator 

- {o"’ 

and the extended spaces 

:= {i; : [0, cx)) ^ K” : PtV G Li VT G [0, oo)}; 

:= {n : [0, oo) ^ M” : PtV G Li+ VT G [0, oo)}. 

In what follows, the superscript n is often suppressed for 
notational simplicity. A (nonlinear) operator A : Lie ^ Lie 
is said to be causal if PtAPt = PtA for all T > 0. A 



causal A is called bounded if the Lipschitz bound 

||PTAa;||i 


l|A|| := 


sup 

T>o-\\PTx\Ui.Q lKTa;||i 

A is said to be monotone [1] if for every mi,M 2 £ Lie, 
Ul > U2 =► Aui > Au2, 


< 00 . 


( 1 ) 


where the inequality is taken coordinate-wise, that is ui > U 2 
means each coordinate of ui is larger than or equal to the 
corresponding coordinate of U 2 - In the case where A is linear, 
monotonicity is equivalent to ALie+ C Lie+. Note that a 
monotone nonlinear operator A which satisfies AO = 0 also 
satisfies ALie+ C Lie+. Such a system is referred to as 
a positive system in this paper. A positive A is said to be 
bounded on Lie_|_ if 


l|A||+:= 


sup 

a;GLie+;T>0;||FTa:|| 17^0 




< 00 


( 2 ) 


The space of bounded linear causal operators mapping 
from Li to itself is denoted by j5f(Li,Li). Each operator 
in j5f(Li,Li) has a natural causal extension to an operator 
mapping from Lie to Lie, where its Li-to-Li induced norm 
is equal to the Lipschitz bound of the extension [9, Section 
2.4]. Note that G G -Sf(Li, Li) is positive if, and only if, its 
causal extension satisfies GLie+ C Lie+. 


Ul 

Gi 









U 2 


►0 


Fig. 1. Positive feedback interconnection of positive systems. 


Consider the positive feedback interconnection of nonlin¬ 
ear causal Gi : Lie —> Lie and G 2 : Lie ^ Lie illustrated 
in Figure 1. Mathematically, 

f U2 = GiMi -|- d2 
^ Ml = G 2 M 2 -l- di- 

Denote the feedback interconnection of Gi and G 2 by 

[Gi, G2], 

Definition 2.1: [Gi,G 2 ] is said to be well-posed if the 
map (mi,M 2 ) H- {di,d 2 ) defined by (3) has a causal inverse 
H on Lie. It is positive it is well-posed and the inverse H 
is positive. It is stable on Lie+ if it is well-posed, positive, 
and the inverse H is bounded on Lie+ (cf. (2)). It is stable 
if it is well-posed and the inverse H is bounded (cf. (1)). 

Note that if Gi and G 2 are bounded and G 2 is linear as is 
the case for standard IQC analysis [ 8 ], stability of [Gi,G 2 ] 
is equivalent to (7 — G 2 Gi)“^ being bounded and causal. 

Remark 2.2: The notion of stability on Lie+ is introduced 
here to characterise the type of well-posed feedback inter¬ 
connections of positive systems (cf. Figure 1) that remain 
stable under disturbances di,d 2 C Lie+. In other words. 


feedback stability is achieved as far as positive signals are 
concerned. The next lemma shows that for linear systems, 
this is equivalent to stability on Lie. 

Lemma 2.3: A linear causal positive A : Lie —> Lie is 
bounded on Lie+ if, and only if, it is bounded on Lie. 

Proof: Sufficiency is trivial. For necessity, note that any 
X G Lie can be written as a; = X 1 —X 2 , where a:i, 0:2 G Lie+ 
satisfy IIPt^iHi < ||PTa;||i and ||Pra: 2 ||i < \\Ptx\\i for all 
T > 0. Moreover, since A is linear. Ax = Aa:i — Aa; 2 , 
whereby 

IIPtAxIIi = ||PTAPTa:||i 

< ||PTAPTa:;i||i -I- ||PTAPTa; 2 ||i 
<||A|| + ||PTXi||i + ||A|| + ||Pra; 2 ||i 
<2||A|| + ||Pra:||i. 

This establishes the boundedness of A on Lie. ■ 

B. Feedback interconnections of positive LTI systems 

Given a matrix A G G denotes its 

transpose. The inequality A > 0 (A > 0) means that 
all elements of the matrix X are positive (nonnegative). A 
matrix A G is called Hurwitz if all its eigenvalues 

have negative real parts. It is said to be Metzler if its off- 
diagonal entries are all nonnegative, i.e. Aij > 0 for i 7 ^ j. 
The spectmm of a square matrix M G and spectral 

radius are defined, respectively, by 

spec(M) := {A G C : A7 — M is a singular matrix}; 
p{M) := max{|A| ; A G spec(M)}. 

Proposition 2.4 ([5, Prop. 1 & 2]): Given a Metzler A G 
jjraxra, ijjg following statements are equivalent: 

(1) A is Hurwitz; 

(ii) There exists a 2 ; G M” such that 2 ; > 0 and z'^A < 0; 

(hi) —A~^ exists and has nonnegative entries. 

Given a P G the following statements are equivalent: 

(iv) p{B) < 1; 

(v) There exists a z G K” such that z > 0 and z^B < z'^; 

(vi) (7 — B)~^ exists and has nonnegative entries. 

Some preliminary results on the feedback interconnection 

of two LTI systems are provided below. They demonstrate 
the fact that static-gain conditions play a crucial role in the 
robust feedback stability of positive LTI systems. Suppose 
that Gi, G 2 have real-rational proper realisations: 

Gi(s) = Gi(s7-Ai)-iPi+Pi and 
G2{s) = C2{sl - A2)-^B2 + P 2 , 

where Ai and A 2 are Metzler and Pi > 0, P 2 > 0, Gi > 0, 
G 2 > 0, Pi > 0, and P 2 > 0. These imply that Gi and G 2 
are positive [5]. In the following the " notation will be used 
to denote the equivalent frequency-domain representation of 
an LTI system or the Laplace transform of a signal in Li as 
above. 

Proposition 2.5 ([10, Thm. 4]): Given Gi and G 2 with 
realisations as in (4), if P 1 P 2 —7 (or equivalently, P 2 P 1 —7) 
is Hurwtiz, then [Gi,G 2 ] is positive. Under this condition. 




[Gi, G 2 ] is stable if, and only if, Ai, A 2 , and Gi( 0 )G 2 ( 0 )—/ 
(or equivalently, G 2 ( 0 )Gi( 0 ) — I) are all Hurwitz. 

Corollary 2.6: Given G\ and G 2 with realisations as in 
(4), if p{DiD 2 ) < 1, then [Gi,G 2 ] is positive. Under this 
condition, [Gi,G 2 ] is stable if, and only if, Ai and A 2 are 
Hurwitz and p(Gi( 0 )G 2 ( 0 )) < 1. 

Proof: First note that repeated applications of Proposi¬ 
tion 2.4 yields that the Hurwitzness of D 1 D 2 —I is equivalent 
to [I — DiD 2 )~^ > 0 , which in turn is equivalent to 
p{DiD 2 ) < 1. Also by Proposition 2.4(iii), Gi( 0 )G 2 ( 0 ) > 0 
whenever Ai and A 2 are Hurwitz. By the same token as 
before, the Hurwitzness of Gi( 0 )G 2 ( 0 ) — / is equivalent to 
/ 5 (Gi( 0 )G 2 ( 0 )) < 1. The claim thus follows from Proposi¬ 
tion 2.5. ■ 

Denote by Hqo the Hardy space of holomorphic trans¬ 
fer functions in the right-half complex plane, which admit 
equivalent LTI operators mapping from the space of square- 
integrable functions 

L2 := I?; : [ 0 ,00) —> 00 : It'COb dt < 00 

to itself via the Laplace transform isomorphism [11, Thm. 
3.30], where [ • [2 denotes the Euclidean norm. 

Proposition 2.7 ([12, Thm. 3]): The feedback intercon¬ 
nection of positive Gi, G 2 G Hqo is such that (/ — 
G 2 Gi)~^ G Hoc and (7 — G 2 Gi)“^ is positive on L 2 , if 
and only if p(Gi( 0 )G 2 ( 0 )) < 1 . 

Remark 2.8: While Proposition 2.7 is different from the 
setting of this paper in that the former considers the L 2 space 
as opposed to Li, it serves to illustrate the important fact that 
the robust stability of a positive closed-loop is determined by 
the static gains of the open-loop systems. Proposition 2.7 can 
also be seen as a generalisation of aspects of Corollary 2.6 
to distributed-parameter transfer functions in Hqo, which do 
not admit realisations of the form G(s7 — A)~^B + D. 

In Section IV, it will be demonstrated that the condition 
/ 5 (Gi( 0 )G 2 ( 0 )) < 1 can be generalised via integral linear 
constraints for guaranteeing robust closed-loop stability of 
positive systems that are time-varying or nonlinear. 


III. Positivity of nonlinear feedback systems 


This section demonstrates that a positive feedback inter¬ 
connection of two positive causal nonlinear (not necessarily 
bounded) systems can be guaranteed to be positive by the 
same condition by Willems [9] that ensures the feedback is 
well-posed. Consider the feedback interconnection of Gi : 
Lie —> Lie and G 2 : Lie —> Lie in Figure 1. It is assumed 
that for i = 1,2, 

(i) Gi is causal and positive; 

(ii) Gi is locally Lipschitz continuous on Lie, i-C- 


sup 

x.ye'LiclPTXjtPTV 


WPAGjX 

WPrix 


Giy)h 

y)\\i 


< 00, VT > 0 . 


The instantaneous gain of Gi is defined by 


a{Gi) 


sup inf sup 

T>0 ‘AT>0 = 

-^T + AT(®~y)7^0 


\\PT+AT{GiX 

\\Pt+at{x 


G^y)h 

y)\\i 


Note that a{G) = 0 if G is strictly causal, i.e. the input 
has no instantaneous effect on the output. It is shown in [9, 
Thm. 4.1] that if a(Gi)a(G 2 ) < 1, then the [Gi, G 2 ] is well- 
posed as per Definition 2.1. The following demonstrates that 
a(Gi)Q;(G 2 ) < 1 also guarantees positivity of [Gi,G 2 ]. 

Theorem 3.1: If a(Gi)a(G 2 ) < 1, then [Gi,G 2 ] is well- 
posed and positive. 

Proof: Well-posedness is established in [9, Thm. 4.1]. 
As for positivity, let u := {ui,U 2 ), d := {di,d 2 ), and 
G(ui,U 2 ) := (G 2 U 2 ,GiUi). The feedback equation (3) can 
be written as 

u = Gu + d. 

Well-posedness and positivity of [Gi,G 2 ] is thus equivalent 
to the causal positive invertibility of (7 — G). [9, Thm. 4.1] 
establishes that when a(Gi)a(G 2 ) < 1 , Pto+atGPto+at 
is contractive for Tq > 0 and sufficiently small AT > 0, 
whereby Pto+at{I — G)Pto+at is invertible. Invertibility 
of 7 — G is then shown on consecutive intervals that cover 
the whole half line [0, 00 ). Moreover, causality of (7 — G)“^ 
is established from the fact that G is causal and the invert¬ 
ibility of 7 — G is based on the convergence of successive 
approximations 

Un+l = GUn + d n = 0,1,... 

for arbitrary uq. Therefore, when G is positive, d > 0, 
and Uq > 0, it follows that lim Un > 0. This establishes 

n—foo 

positivity of (7 — G)“^ along the same lines of causality. ■ 
Consider again the setting of Proposition 2.5, where 
Gi(s) = Gi(s7 - Ai)-^Bi + Di and G 2 (s) = G 2 (s 7 - 
A 2 )~^B 2 +D 2 with Ai and A 2 Metzler and Bi > 0, B 2 > 0, 
Gi > 0, G 2 > 0, Di > 0, and D 2 > 0. Observe that 

a{Gi) = sup \Diu\ =: ||A||i^i- 

|u| = l 

Thus, Theorem 3.1 ensures that [Gi, G 2 ] is positive whenever 
|| 79 i||i^i|| 7 ) 2 ||i^i < 1. Alternatively, it may be established 
by the first part of Corollary 2.6. To see this, note that 

p(79i7?2) < ||7)i7)2||i^i < ||7?i||i^i||772||i^i < 1- 

Theorem 3.1 is useful in that it can be applied to nonlinear 
time-varying systems. 


IV. Integral linear constraints 


Given a causal bounded A : Lie —> Lie, define the graph 
of A with respect to positive signals as 



G Li+ :y = Ax 


Y 


Similarly, define the inverse graph as 



G Li+ : X = Ay 


Y 


The following definition of integral linear constraint follows 
that in [3]. As will be seen later, this allows the establishment 
of a link from the main robust stability result to the static 
gain condition in Corollary 2.6. 



Definition 4.1: A causal bounded system A : L™ —> 
is said to satisfy the integral linear constraint (ILC) defined 
by the multiplier H G Mix(m+p) 

nw(f) dt > 0 Vw G 0(A). 

This is denoted A G ILC(n). On the contrary, A is said to 
satisfy the strict complementary ILC if 

pOO pOO 

/ nu;(f) dt <—e / dt Vw G 0'(A) 

■Jo Jo 

for some e > 0, where G K” denotes the vector whose 
entries are all ones. This is denoted A G ILC‘^(n). 



Remark 4.2: Having 11 as a bounded LTI operator from 
Li to Li does not generalise the definition of an ILC for 
the following reason. Denote by tt the convolution kernel 
corresponding to 11 and * the convolution operation, then 



= 7r(0) 


■ m' 
_^(0) 


(f) dt 
dt 


s=0 


That is, only the static gain of tt, which is a constant matrix, 
matters in the value of the integral. 

Lemma 4.3: Suppose A : —> L^’g is a causal bounded 

positive LTI system, then A G ILC(n) is equivalent to 


n 



> 0 , 


(5) 


where A denotes the transfer function representation of A. 
Similarly, A G ILC‘^(n) is equivalent to 



n 

A(0) 

Proof: For any u G Li+, 

i 

•00 

n 

1 

u{t) 

(Au)(t) 

pOO 

= / n 

.70 L 

u(t) 

(Au)(t) 

= n 

A(o). 

u(0)- 

note that 11 

A(o) 

u(0) 


< 0 . 


dt 


e dt 


s=0 


m( 0) > 0 for all u G Lm_ is 


equivalent to 11 


4 
A(0). 


> 0, since u(0) = u{t) dt > 0. 

The second part of the claim can be shown using the same 
lines of arguments. ■ 


Note that the static-gain condition (5) in Lemma 4.3 is 
amenable to distributed verifications, since it only involves 
multiple vector multiplications. The computation also scales 
linearly with the dimensions of the system concerned, much 
in the spirit of [5]. 


V. Robustness analysis 


The main robust stability result for positive feedback 
systems is stated. Note that the open-loop systems do not 
need to be positive for the following to hold. 

Theorem 5.1: Given a (nonlinear) bounded causal Gi : 
—> Lj’g and a linear bounded causal G 2 : L^’g ^ 
suppose there exists a 11 G such that 

(i) [tGi, G 2 ] is well-posed and positive for all t G [0,1]; 
(11) tGi G ILC(n) for all r G [0,1]; and 
(ill) G2GlLC'=(n). 

Then [Gi,G 2 ] is stable on Lie+. In addition, if Gi is also 
linear, then [Gi, G 2 ] is stable. 

Proof: By hypothesis, for all r G [0,1], v G Q{tGi) 
and w G G'{G 2 ), there exists e > 0 such that 



nri(f) dt > 0 


and 


pOO pOO 

/ nu;(f) dt < —e / lp_f_^w{t) dt. 

Jo Jo 


Define T* := 211-1-The inequalities can be written as 

pOO pOO 

/ ^v{t)dt>e = e||n||i 

Jo Jo 


and 


pOO pOO 

/ ^w{t)dt<-e l'^_^_^w{t) dt =-e\\w\\i, 

Jo Jo 

where the fact that v, w G Li+ has been exploited. Therefore, 
by the feedback configuration in (3), 

/ OO 

^{v{t) - w{t)) dt 


< 




- 01 , 


dt 


= V'dMi 


di (f) 

-d2{t) 

T di{t) 

P+^ [d2{t) 

I|rf2||l),' 


dt 


for all V G 0(tGi), w G G'{G 2 ), and di,d 2 G Li+, 
where := max Ifiil and T' = [0i,... fim+p]- In other 

words. 


IHIi + lhlli<^(||di||i + ||rf2 


( 6 ) 


Note that when d 2 = 0, di = (/ — tG 2 Gi)ui by linearity 
of G 2 . As such, it follows from (6) and v = 
for all ui e Li-h, 


Ui 

rGiUi 


that 


lluill < |||(/-rG2Gi)ui||i. (7) 

By the well-posedness and positivity assumption, the in¬ 
verse (/ — tG 2 Gi)~^ is well-defined on Lie and satisfies 
(/ — rG 2 Gi)“^Lie+ C Lie+. Given any u G Lie+ and 
i' G [0,1], define 


UT := (/ - i^G2Gi)-^Pt{I - i^G2G0u G Lie+. 



Then 

||^tu||i = ||PtWt||i 

^ ~ 1^G2Gi)ut\\i 

= ^||Pt(/ - rG 2 Gi)u + (r - v)PtG 2 Giu\\i 

< ^||PT(/-rG2Gi)u||i 

+ ||r-HI|G2||||Gi||||PT«||i, 

( 8 ) 

where (7) has been used to arrive at the second inequality. 
Since (7 — i^G 2 Gi)“^ is bounded on Lie+ for u = 0, it 
follows from (8) that (7 — tG 2 Gi)“^ is bounded on Lie+ 
for all T < (^||G 2 ||||Gi||)“^. By an inductive argument and 
repeated applications of (8), it follows that (7 —rG 2 Gi)“^ is 
bounded on Lie+ for all t G [0,1]. In other words, [rGi, G 2 ] 
is stable on Lie+ for all t G [0,1], 

For the final claim, note that linearity of Gi and G 2 
implies that of (7 — G 2 Gi)“^ and apply Lemma 2.3. ■ 

VI. Examples 

A. LTl systems 

Consider the feedback interconnection of two positive LTI 
systems Gi and G 2 , where Gi(s) = Gi(s7 — Ai)~^Bi + 
Di and G 2 (s) = C 2 {sl — A 2 )~^B 2 + D 2 with Ai and A 2 
Metzler and Bi > 0, B 2 > 0, Gi > 0, G 2 > 0, Di > 0, and 
D 2 > 0. Suppose that ||7?i||i^i||752||i^i < 1, whereby 

||r77i||i^i||772||i^i < IVtG [0,1], 

so that [rGi, G 2 ] is well-posed and positive by Theorem 3.1. 
It is shown below that the sufficiency part of Corollary 2.6, 
i.e. p(Gi( 0 )G 2 ( 0 )) < 1 implies that [Gi,G 2 ] is stable, can 
be recovered from the main result Theorem 5.1 using a 
particular multiplier 11. 

To begin with, note that p(Gi(0)G2(0)) < 1 is equivalent 
to the existence of a z > 0 such that 2 ;^(Gi( 0 )G 2 ( 0 )—7) < 0 
by Proposition 2.4. Define 

n:=z^[Gi(0) -7]. 

It is straightforward to see that 

n *^2^0) I >0Vtg[0,1]. 

1 J [rGi(UjJ 

By Lemma 4.3, these are equivalent to rGi G ILC(n) for 
all r G [0,1] and G 2 G ILC‘^(n). As such, [Gi, G 2 ] is stable 
by Theorem 5.1. 

B. Feedback channels with time-varying gain 

Suppose two strictly causal bounded positive LTI systems 
Ml and M 2 are connected in a positive feedback loop via 
channels with time-varying gains 5il and 82 ! satisfying 0 < 


Si{t) < a and 0 < (52(f) < P- Stability of this is equivalent 
to the stability of the stmctured feedback interconnection 
[M, A], where 



Ml' 

0 


and A = 


■<5i7 

0 



Note that with 11 := (a7,/37, —7, —7), tA G ILC(n)VT G 
[0,1]. By the strict causality assumption on M, application of 
Theorem 3.1 yields that [M, tA] is well-posed and positive 
for all T G [0,1]. It thus follows from Lemma 4.3 and 
Theorem 5.1 that [M, A] is stable if 



VII. Conclusions 

The positivity of positive feedback interconnections of 
positive nonlinear systems is shown to be guaranteed by 
an instantaneous gain condition. The static gain conditions 
on robust feedback stability of positive LTI systems have 
been generalised to time-varying and nonlinear systems via 
the notion of integral linear constraints. Future research 
directions may involve investigating if the use of linear 
time-varying multipliers reduces conservatism in robustness 
analysis and the use of Loo based conditions as a replacement 
of the Li based ILC. 
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